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1.  INTRODUCTION 

The  prediction  of  radar  and  ESN  performance  in  a  complex  mari¬ 
time  environment  hat  become  a  routine  task  for  the  naval  operations 
analyst.  With  the  emphasis  on  operational  rather  than  environmental 
variables,  resort  is  nsnally  had  to  readily  available  computer  programs 
for  propagation  and  detection.  If  these  are  not  available,  several 
physical  and  mathematical  approximations  mnst  be  made  to  facilitate 
routine  predictions. 

Simple  ladar  models  employing  'standard'  environment  and  detector 
have  proved  useful  for  ad  hoc  estimates  of  performance  nnder 
globally  averaged  conditions.  These  models  however  often  employ  siaqile 
approximations  for  multipath  effects,  signal  statistics,  detection 
criteria,  attenuation  and  backscatter  and  care  is  also  needed  to  ensnre 
that  inpnt  variables  are  not  outside  the  range  of  validity  of  the 
model . 


The  goal  of  the  radar  model  described  herein  was  threefold  - 
(i)  it  should  produce  reliable  data  over  a  wide  range  of  environmental 
and  radar  parameters,  excluding  ducting  conditions,  (ii)  the  programs 
should  be  written  in  a  high-level  scientific  language  and  (iii)  the 
model  should  be  programmable  on  medium  sixed  mini  and  microcomputers. 
In  reference  1  a  BASIC  program  was  described  which  was  specifically 
designed  for  the  Tektronix  4051  system,  sad  which  was  concerned  with 
the  performance  of  scanning  search  radars  with  fixed  threshold  detectors. 

This  model  has  since  been  enlarged  and  rewritten  in  PASCAL,  and 
the  detection  subroutines  have  been  extended  and  modified  to  improve 
numerical  stability.  Both  fixed  threshold  and  (TAR  detection  are 
treated  for  coherent  and  non-coherent  detectors  and  various  n-out-of-n 
detection  criteria.  This  paper  also  presents  a  more  detailed  description 
of  the  physical  principles  and  range  of  validity  of  the  modal.  Each  of 
the  major  subroutines  is  covered  by  a  separate  annex.  To  assist  with 
program  development,  the  results  of  eaeh  routine  are  also  presented  in 
graphical  form. 

With  many  aspects  of  radar  performance  prediction  the 
literature  presents  several  alternative  teehai«ues  or  relations  used  hy 
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various  authors  with  qualified  reliability.  The  most  important  of  such 
topics,  which  are  considered  to  merit  separate  annexes,  are 

(i)  refractivity  and  the  effective  earth  radius 

(ii)  attenuation  by  the  atmosphere 

(iii) pattern  propagation  factor 

(iv)  antenna  pattern  function 

(v)  antenna  noise  temperature 

(vi)  surface  and  volume  clutter 

(vii) calculation  of  detection  and  acquisition  probabilities 

The  annex  on  clutter  is  predominantly  concerned  with  a 

description  of  the  surface  clutter  routines,  which  have  not  been 
reported  elsewhere.  Probability  of  detection  algorithms  have  been 
validated  for  accuracy  and  nuawrical  stability  over  a  wide  range  of 
probabilities  and  target  scintillation  parameters. 

2 .  THEORETICAL  BACKGROUND 

The  theory  of  radar  detection  is  well  documented  in  the 

literature  (e.g.  references  2-7),  but  will  be  briefly  reviewed  below 
to  give  a  broad  outline  of  the  method  used  and  to  place  in  context  the 
symbols  and  terms  used  in  later  sections. 

2,1  The  Radar  Equation 

In  free  space  electromagnetic  radiation  spreads  spherically,  so 
that  the  power  density  at  distance  R  in  a  tangential  surface  of  area  A 
is 


P,  P»6» 

X  *  W  (  1  > 

where  Pt  is  the  transmitted  power  and  Gt  is  the  power  gain  of  the 
transmit  antenna  in  the  direction  of  propagation.  Iqnation  1  is  the 
one-way  transmission  equation  appropriate  to  conniuaicat ions  or  E8M/BCM. 
In  the  latter,  A  can  be  replaced  by  the  effeetive  antenna  area  A,  to 
caloulate  the  received  power  in  terms  of  the  wavelength  (X)  and  the 
receive  antenna  gain  (Of) 
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(  2  ) 


The  'radar  equation’  in  its  simplest  form  assumes  that  the 
radiation  incident  on  a  scatterer  is  captured  in  an  area  a  and 
reradiated  isotropically.  This  yields  the  bistatic  radar  equation 


PtGrGtX*c 


(  3  ) 


where  R,  is  the  transmit  antenna  to  target  range,  and  R,  is  the  target 
to  receive  antenna  range.  For  monostatic  radar  Gr=Gt=G  and  R1*R1=R. 
If  more  than  one  path  is  significant  the  received  field  is  modified  by 
a  pattern  propagation  factor  F.  The  usual  definition  of  F  (Annex  C) 
includes  the  antenna  pattern  functions  (Annex  D)  of  the  direct  and 
indirect  rays,  so  that  G  is  replaced  by  the  gain  along  boresight  (G«) 


PtG.*X*F«o 
Pr  -  Tfa  »l*L' 


(  4  ) 


The  propagation  path  is  calculated  using  the  rules  of  geometric 
optics  (ray  theory)  for  a  smoothly  decaying  pressure  and  moisture 
gradient,  with  the  geometric  variables  transformed  to  the  effective 
earth  radios  model  to  simplify  calculations.  The  approach  used  in  this 
model  to  calculate  the  pressure  and  humidity-dependence  of  the 
effective  earth  radius  is  described  in  Annex  A. 


Pattern  propagation  factors  for  targets  over  the  horixon  are 
evaluated  using  diffraction  theory  (ref  3),  and  an  interpolation 
procedure  is  used  for  targets  near  the  horlson. 

Losses  due  to  absorption  by  atnospherlo  gases  and  precipitation 
are  accounted  for  by  the  loss  factor  L.  This  loss  is  computed  using 
the  calculated  attenuation  rate  (Annex  I)  and  the  slant  range.  System 
losses  are  not  included  in  aquation  4  slnee  the  power  is  referred  to 
the  antenna.  The  incorporation  of  additional  losses  due  to  filter 
mismatoh,  pulse  compression,  etc  are  deseribed  in  reference  1. 

2.2  C lit tar 

Thu  clutter  return  from  the  sea  surface  (P#)  is  calculated  from 
the  radar  aquation  (eqm  4)  with  the  olutter  cross-sect  ion  given  by 
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a  =  o»A  =  a*  lfi|§-*ecr  (  5  ) 

where  (P  is  the  effective  horizontal  beamwidtb,  c  is  the  speed  of 
light,  y  1*  the  grazing  angle  and  B  is  the  bandwidth  of  the  matched 
filter.  Expressions  for  the  clntter  cross-section  per  cnit  area  (o*> 
ere  described  in  annex  F. 


Volnae  backscatter  from  rain  and  snow  is  determined  in  an 
analogous  manner,  aith 


o 


avV 


(  6  ) 


2.3  K'oise 


Quasi-thermal  rcise  is  evaluated  from  the  system  noise  temp¬ 
erature  (Tg)  referred  to  the  antenna.  The  equivalent  noise  power 
(Pn)  is  then 

P„  «  E  /x  -  kT./r  (  7  ) 

n  n  s 

where  x  is  the  radiated  pulse  width  and  ^  is  the  spectral  density  of 
the  noise  power.  (Care  should  be  taken  using  bandwidth  in  the  radar  equation 
if  x  is  not  of  order  1/Bn) .  The  system  noise  temperature  is 

T,  «  T,  +  Tt(Lr-l)  +  LrT,(NF-l)  (  8  ) 

in  which  Tg  is  the  antenna  temperature,  Lr  and  Tr  are  the  receiving 
line  losses  and  temperature  respectively,  T("290  K  and  NF  is  the 
zeceiver  noise  figure.  The  method  for  evaluating  T(  is  outlined  in 
annex  E. 


2.4  Probability  of  Detection 

Gaussian  statistics  are  assumed  for  elutter  and  noise,  and  the 
single  pulse  signal-to-noiae  ratio  is  defined  as 


S/N  -  P,/(P«  ♦  P„)  <  9  ) 

t  6  * 

The  probability  of  paint  dependa  on  S/N,  the  probability  of 
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false  alarm,  the  degree  of  pre-detection  (coherent)  and  post-detection 
(non-coherent)  integration,  and  the  target  fluctuation  statistics.  A 
generalized  Chi-square  distribution  is  used  to  synthesize  the  amplitude 
statistics  of  the  target  return,  with  a  single  distribution  paraveter 
which  is  determined  by  the  relative  timescale  of  target  scintillation. 
Non- fluctuating  targets,  Swerling  cases  I-IV  and  Veinatock  targets  ere 
included.  Programs  for  evaluating  paint  probability  from  S/N  are 
described  in  annez  6.  A  general  method  for  deteia<ining  acquisition 
ranges  for  m-ont-of-n  detection  criteria  it  also  described  together 
vith  results  for  1/1,  2/2.  and  2/3  detectors. 
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ANNEX  A 

REFRACT! VITY  AND  EFFECTIVE  EARTH  RADIUS 

Variation  in  the  refractive  index  (n)  of  air  with  altitude 
causes  electromagnetic  rays  to  follow  curved  paths  through  the 
atmosphere.  The  classical  Method  for  dealing  with  this  phenomenon,  when 
calculating  ray  geometry,  is  to  replace  the  earth*  a  radius  (a)  with  an 
effective  earth'a  radius  (ae) 


ae  *=  Ka 


(A-l) 


and  to  assume  an  equivalent  homogeneous  atmosphere.  From  Snell's  Law 

K  =  [1  +  a^-1  (A-2) 


where  gg*  is  the  change  in  refractive  index  with  ^e-ght.  Replacing 
refractive  index  with  refractivity  N  (*  10#(n~l))  gives 


[1 


-i 


(A-3) 


The  model  makes  the  further  assumption  that  N  varies  exponentially  from 
its  sea  level  value  Ng  to  its  value  at  the  tropopanae  Nt 

V°«etf  (A-4) 

N  <=  N  e  *  * 


and 


Na 

Tf  lo*. 


ft*"** 


(A-5) 


where  h{  is  the  height  of  the  tropopanae.  Since  moat  of  the  ray  bending 
oconra  in  the  firat  few  kilometrea  of  the  atmoaphere,  where  '££■  ia 
relatively  conatant,  a  linear  approximation  it  made  naing  the  value  of 
tea  level. 


r  •  N.  N.1-1 

y-  *  i0*hjlo«eiq 


(A-6) 


and  the  tropopanae  ia  atsemed  to  have  a  conatant  height  of  IS  kilometers 
with  the  refractivity  conatant  at  61  uaita  at  that  height.  Snbatltnting 
for  Nt  and  ht 


.../A-2 


K  =  ^1  +  0.00049  N$  logg-^-j'1 


The  refractive  index  of  a  pore  gas  is  given  by  the  Debye  relation 
(ref  8)  and,  assuming  ideal  gas  behavionr,  this  becomes 


where  P  is  the  pressure,  R  is  the  gas  constant,  a  is  the  electronic 
polarizability,  p  is  the  electric  dipole  moment,  T  is  the  temperature  and 
ee  is  the  permittivity  of  free  space.  For  a  mixtnre  of  gases  P  is 
replaced  by  the  partial  pressures  and  the  term  on  the  right  is  summed 
over  all  species  which  contribute  to  the  polarization  -  mainly  oxygen, 
nitrogen  and  water  for  air  near  sea  level.  This  equation  is  valid  only  if 
the  electric  dipole  p  can  rotate  in  phase  with  the  applied  electric  field. 
At  RF  and  low  microwave  frequencies  and  for  a  normal  atmosphere  this 
reduces  to 


,,  4810P-  1 

N  -  10 ‘  (n-1)  -  Z^|p  +  - 5:-*  I 


where  P  and  Pw  are  total  pressure  and  partial  pressure  due  to  water 
vapour  in  mbar  and  T  is  the  temperature  in  Kelvin 


The  lowest  frequency  rotational  resonance  for  water  is  at  22  GHz 
( f 0) .  Using  the  Krarers  -  Konig  relationship  and  Lorentxian  line  shape 
sssumption  (i.e.  collisions!  broadening),  the  variation  in  the  real  part 
cf  the  dielectric  constant  across  an  absorption  band  is  given  by 


6e*  *  1m  lo.ifc 


(A- 10) 


where  y***  i»  the  peak  absorption  (dB/km),  «  is  the  frequency  and  c  is 
the  speed  of  light.  Taking  a  worst  ease  (100%  humidity  and  300  K)  8»s  “ 
0.4  x  10~6  (or  AN  -  0.2  for  N  «  400)  across  the  22  GBz  absorption  band  - 
an  insignificant  error.  The  general  variation  in  N  due  to  disperaion  in 
the  molecular  polarixability  from  the  IR  to  RF  regiion  (ref  8)  could 
account  for  10  -  40%  in  N.  However,  since  the  range  of  validity  of  the 
overall  model  is  restricted  to  RF  and  mlorowave  bands,  only  a  mean  value 
for  refractivity  dispersion  is  used  is  the  model. 

The  model  cses  environmental  inputs  of  temperature,  total 


.../A-S 


pressure  sod  relative  humidity  (H)  to  calculate  Ns  (using  equation 
(A-9)).  To  calculate  partial  pressure  due  to  water  vapour  (fro*  the 
input  parameters)  an  eapirical  {ora  of  the  Clausius  -  Clapeyron  relation 
is  used 


(A-U) 


where  kj  =  1.8178  i  10’  abar 
k,  »  5329  l 

and  H  is  the  relative  humidity  in  percent 


If  the  relative  humidity  is  not  known  explicitly,  equation  (A~ll) 
can  be  rewritten  in  terms  of  wet  and  dry  bnlb  temperatures.  The  saturated 
vapour  pressure  at  the  wet  bulb  temperature  ia  obtained  from  equation 
(A-ll)  with  T  *  Tw  and  H  *  100.  To  obtain  the  vapour  pxeaaure 
at  aabient  teaperaure,  the  vet  bulb  depression  AP,  is  subtracted.  For 
total  presaure  P  in  abar,  T  in  degrees  Celsius,  this  is 

AP„  =  [0.64 6  x  10-*  P  (1  +  0.00094TW)J(T  -  T.)  (A- 12) 

The  considerable  variability  of  N,  and  hence  K  ia  illustrated  in 
figures  A-l  and  A-2.  In  figure  A-l,  N(  ia  plotted  against  H  for  several 
teaperaturea  at  a  constant  pressure  of  1019. 25  abars.  In  figure  A-2,  K  is 
plotted  against  H  for  a  similar  range  of  values. 


AbNEX  B 


AITENPATION  BY  THE  ATHOSPHEttE 

Atmospheric  Absorption  Dae  to  Bncondcnsed  Gases 

Genera 1 


The  one  way  absorption  (Aj),  due  to  constituent  gas  j,  for  a  beam 
traversing  a  path  in  the  ataosphere  is  obtained  by  integrating  the 
absorption  coefficient  (yj)  over  the  ray  path 


S  =  j* Yj  <*’  d* 


(B-l) 


where  a  «  the  distance  along  the  ray  path 

R  *=  the  total  length  ef  the  raypath  (ie  the  slant  range) 


Assuming  absorption  is  proportional  to  partial  pressure  with  an 
exponential  model  for  each  atmospheric  constituent,  the  absorption 
coefficient,  as  a  function  of  height,  is 

-h/H, 

T j (h)  -  Tj0  e  i  (B-2) 

where  yj0  is  the  coefficient  of  absorption  at  sea  level  fcr  constituent  j 
and  Bj  is  a  constant  for  constituent  j.  The  height  of  the  direct  ray  at 
a  distance  (a)  along  the  ray  path  is  given  by 

h  »  J s*+2s(a#+h»)ain#e+(a#+hj)*'  -a,  (B-3) 

where  a,  ■  effective  earth's  radius 

ht  -  antenna  height  above  sea  level 

I,  •  elevation  angle  of  the  ray  at  the  antenna 


Subatituting  equation  B-2  into  equation  B-3  gives 

/s1'»^s(a>eh>)si«#,'»(a,'»hi)»t  -a, 

tj(.)  -  Tj0« 


<B-4) 


>nd  therefore 


VS)  =  Yj°fe 


/*1+2*(«.+hi)  sin#_+(te+hi)  *  -a. 


This  integral  is  solved  numerically  in  the  model  nsing  n-point 
Gaussian  quadrature  techniques  in  which 


Aj  =  yjoR  f(xi'R) 


Yt»+2sUe+hi)  »ln#e+(ie+hi)*  -ag 
where  f(s)  «  e  j 


and  Xj  and  Wj  are  tabulated  in  column  k=0  of  table  25.8  of  reference  15 


for  various  n.  In  this  model  n*=8  is  used. 


Total  Absorption 


The  total  one  way  absorption  (A)  is  calculated  by  summing  Aj  for 
all  constituent  gases  where  A  and  Aj  are  in  decibels: 

A  -  ^Aj  (B-8) 

j 

At  B'icrowave  frequencies,  absorption  of  electromagnetic  radiation 
by  nncondensed  gases  is  due  to  the  interaction  with  sioleonlar  dipoles 
rotating  on  a  10*  Ht  timescale  (ref  9).  The  strongest  molecular  absorber 
is  water  which  has  a  permanent  electric  moment  and  gives  rise  to  a 
strong  line  at  22.235  GHz  and  other  lines  in  the  millimetre  region. 
Oxygen  lacks  s  permanent  electric  moment  bnt  displays  a  magnetic  fine 
structure  at  120  GHz  and  some  fifty  observable  lines  between  40  and  80 
GHz.  Although  this  magnetic  interaction  is  much  weaker  than  the  electric 
dipole  interaction,  oxygen  (by  virtue  of  its  arach  greater  denaity)  is 
responsible  for  an  absorption  cosqparable  with  water  at  sea  level  and 
greater  at  higher  altitudea.  No  other  polar  or  paramagnetic  tropospheric 
gases  give  rise  to  significant  absorption  at  the  frequencies  of  interest 
and,  therefore,  in  this  model  total  one  way  abaorption  for  uneondenaed 
gases  is  assumed  to  be 

A  ■  i#  ♦  1,  (B-9) 

A0  and  Aw  are  calculated  from  aquation  6  above  using 


.../B-S 
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B0  =  6.1986  1 i lone t res 
=  2.4261  1  lionet res 

and  values  for  y00  and  ywo  calculated  fro*  the  models  described  below 
using  set  level  environmental  parameters  at  input. 

Absorption  Coefficients  -  General 

Microwave  spectral  'lines'  have  vanishingly  small  width  only  in 
the  limit  of  zero  pressure.  The  width  of  the  spectral  lines  varies 
inversly  with  the  excited  state  lifetimes  or  dephasicg  times,  which  are 
of  the  same  order  as  the  rate  of  'effective'  molecular  collisions.  At 
sea  level  these  occnr  on  a  10*  Bz  time-scale  for  both  water  and  oxygen 
(ref  10).  Lew  altitude  radars  usually  do  not  operate  near  the  peak  of  an 
absorption  line,  however  the  calculated  residual  intensity  in  the  ’tail’ 
of  the  absorption  is  critically  dependent  on  a  model  of  the  relaxation 
process.  An  exponential  decay  is  assumed,  resulting  in  a  Lorentzlan 
power  spectrum. 

Absorption  Coefficient  -  Tatar 

The  main  water  baud  of  interest  in  the  microwave  region  is  the 
22.235  GHz  line  although  there  is  some  background  absorption  from  the 
180  0Hz  and  320  GHz  bands.  From  ref  11,  the  expression  for  the  line  width 
(Af)  is 

Af  »  1.349x10**  *  0. 20846  <p-*w>  W°"}»  s  (B-10) 

where  P  and  Pf  ere  the  total  pressure  and  partial  pressure  due  to  water 
vapour  in  millibars,  T  is  in  degrees  Kelvin.  Although  water  eollisions 
are  more  'effective'  than  those  of  other  molecules,  the  line-width  is 
approximately  proportional  to  the  total  pressure.  The  expression  for 
water  vapour  absorption  is  then  given  by 

Yw  -  1.9014x10'*  jfPw(^]T/* a2* 144 ♦  rrM  dB/kn  <»-ll> 

where  f  le  the  frequency  in  gigaherts  and  F  is  a  shape  factor  given  by 

F  “  ♦  rmffiiiTp]  <*-«> 
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where  f0  =  22.235  GHz  and  Tre«  it  the  residual  vatei  vapour  absorption 
dne  to  the  tails  of  lines  above  100  GHz  given  by  the  empirical  expression 


rres  =  1,591  T-’/*P.Pwf*  dB/km 


(B-13) 


Absorption  Coefficient  -  Oxyxen 

Although  Oj  does  not  have  a  permanent  electric  moment,  it  is 
paramagnetic  and  its  microwave  absorption  is  dne  to  the  fine-structure 
transitions  in  which  the  magnetic  moment  assumes  various  directions  with 
respect  to  the  rotational  angular  momentum  of  the  molecule.  The  coupling 
of  the  rotational  angular  momentum  (quantum  number  N)  with  spin  (quantum 
number  S=l)  results  in  total  angular  momentum  with  quantum  number  J  ■= 
(N-l),  N  or  hi+1.  The  selection  rules  for  the  microwave  trnrsitions  are 


J  “  N  — >  J  «  N  +  1 


(B-14) 


The  relative  intensities  of  these  lines  are  related  to  the  Boltzmann 


population  factors  for  the  J-states. 


Nj  -  (2J+1){|  e 


(B-15) 


Spectroscopic  data  for  0»  (ref  9)  gives  (kT/hB)«145  at  3001. 
To  account  for  contributions  from  overlapping  lines  sumsMtion  up  to 
frequencies  corresponding  to  J«45  have  been  carried  out  to  determine 
the  0a  absorption.  The  positions  of  the  0»  lines  (ref  2)  are  listed  in 
table  B-l.  To  compute  the  total  oxygen  absorption  the  Van  Vleck-Weisskopf 
theory  (ref  12)  is  used  which,  after  some  substitution,  reduces  to 


y0  -  2.005*  PT-'f^jT]*,,  dB/km 
where  AN  -  [F^  ♦  +  F^*].  -2.0f$N(H*l)/T 

Fo  “  p+fnr* 


rN+  "  ITnVTT»*( 


(B-16) 


(B-17) 


(B-lt) 


(B-19) 


(B-20) 


...ns 


LN+l-UffL-ll 


(B-21) 


»*N- 


•*No 


2^mri} 


(B-22) 


•re  line  shape  factors,  P  it  baroaetric  presssre  in  nillibars,  T  is  in 
Kelvin,  f  is  the  frequency  in  GHz  an  d  Af  is  the  line  width. 


To  account  for  the  breakdown  of  the  Van  Vleck  -  Veisskopf  theory 
for  overlapping  lines,  seai-eapirical  line  widths  have  been  need  to  give 
good  agreeaent  at  sea  level.  At  aoderate  altitudes  the  linewidth  narrows 
according  to  the  expression 

Af  =  0.18949  CBz  (B-23) 

Above  SO  ka  doppler  broadening  and  Zeeaan  splitting  aust  alto  be 
inclnded  since  these  begin  to  doainate  the  apparent  linewidth. 


Condensed  Water 


The  attennation  beyond  the  Rayleigh  region  cannot  be  treated  as 
absorption  but  rather  as  a  resonant  anltiple  scattering,  which  depends  on 
the  distribution  of  drop  sizes.  At  very  high  water  densities  such  as  fog 
or  rain,  the  absorption  of  water  can  be  calculated  froa  staple  eapirical 
expressions.  Rain  attenuation  is  usually  expressed  in  the  fora 


Train 


(B-24) 


where  t  is  the  rainfall  rate  in  an/ hr  and  a, fi  axe  functions  of  frequency. 
Rivers  (ref  13)  has  suggested  the  following  seai-eapirical  expreseions 
based  on  Hie  scattering  theory  and  raindrop  aixe  distributions. 

•  “  $l+?^i5f3{MU+l*/i300)  VMl+f*/ 12100) ‘/»  <B_25) 


»  -  1.30  ♦  0.0372  jl-  [l+ 

i.e.  there  is  a  aaxinua  in  the  power  law  at  I -band. 


(B-26) 


Goldstein  (raf  14)  has  deduced  a  foraula  for  attenuation  by  fog: 


) 
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rfog  -  4.89xlO-«Mf*  (B-27) 

where  V  is  water  content  in  g/n».'  Conparing  (B-24)  and  (B-27)  it  is 
clear  that  at  CHF  frequencies  (say  1  GHz)  a  loo/hr  rainfall  rate  gives  an 
attenuation  cooparable  with  a  fog  of  density  0.67  g/o*.  (100  o 

visibility).  Also,  from  eqn  (B-ll).  uncondensed  water  of  the  saoe 
average  density  results  in  about  one  tenth  of  the  attenuation.  Fog 
attenuation  is  often  recorded  as  optical  visibility  (d).  An  eopirical 
relation  between  I?  and  d  can  be  deduced  fron  Goldstein's  data 

K  =  k  d-1-43  (B-28) 


where  k  =  1670  when  d  is  in  feet  and  305e4  when  d  is  in  Metres. 


Table  B-l 

Oxygen  Resonance  Frequencies 


N 

fN+(GHz) 

f  (GHz) 

1 

56.2648 

118.7505 

3 

58.4466 

62.4863 

S 

59.5910 

60.3061 

7 

60.4348 

59.1642 

9 

61.1506 

58.3239 

11 

61.8002 

57.6125 

13 

62.4112 

56.9682 

15 

62.9980 

56.3634 

17 

63.5685 

55.7839 

19 

64.1272 

55.2214 

21 

64.6779 

54.6728 

23 

65.2240 

54.1294 

25 

65.7626 

33.5960 

27 

66.2978 

53.0695 

29 

66.8313 

52.5458 

31 

67.3627 

52.0259 

33 

67.8923 

51.5091 

35 

68.4203 

50.9949 

37 

68.9478 

50.4830 

39 

69.4741 

49.9730 

41 

70.0000 

49.4648 

43 

70.5294 

48.9582 

43 

71.0497 

48.4530 

i 
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ANNEX  C 

PATTERN  PROPAGATION  FACTOR  (P) 


General 


The  form  of  the  radar  equation  used  to  define  the  pattern 
propagation  factor  <F)  is 

P,G*oA»F« 

pr  =  -tenr  (c-i) 

so  that  F  includes  terms  in  the  antenna  pattern  function  and  the  received 
field  doe  to  multipath,  diffraction  and  scattering  effects. 

The  Interference  Region 

In  the  interference  region  (  between  the  radar  and  a  critical 
ground  range  (G0)  ) ,  F  is  calculated  using  ray  theory  for  one  direct  ray 
and  one  ray  reflected  from  the  sea  surface.  Gc  is  near  the  optical 
horizon  for  millimetre  to  centimetre  wavelengths  but  is  typically  only 
60%  of  the  horizon  range  for  UHF  radars.  Expressions  for  the  calculation 
of  Gc  are  at  equations  C-32  and  C-33. 

Within  the  interference  region  (  0<G<Gc  )  F  is  obtained  simply  ss 

F  «  f(6l)yi+x*+2x.cos(5»i)  (C-2) 

where  x  »  (C-3) 

and  f(0)  is  the  antenna  pattern  function  (  see  Annex  D  ) 

8i  is  the  elevation  angle  of  the  direct  ray 
9t  is  the  elevation  angle  of  the  indirect  ray 
9 i  1*  tka  phase  difference  between  the  direct  and  indirect  ray 
r  is  the  roughness  coefficient 
D  it  the  divergence  factor 

p  is  the  magnitude  of  the  sea  surface  reflectivity 

The  phase  difference  (  ^  )  between  the  direct  and  indirect  ray 
is  the  sum  of  the  geometric  phase  difference  due  to  different  path 


./C-2 


C-2 


J 


lengths,  and  the  phase  shift  at  reflection  for  the  indirect  ray  due  to 
dielectric  effects.  Direct  subtraction  of  the  direct  and  indirect  path 
lengths  nay  lead  to  numerical  instability  and  the  preferred  method 
of  calculating  the  path  difference  (6)  is 

£  -  l(#/2n)  •=  (C-4) 

where  the  symbols  are  as  shown  in  figure  C~l.  The  slant  ranges  (Ra,Ra,R) 
are  first  evaluated  as  the  corresponding  ground  ranges  (Ga.Ga.G)  which  are 
the  solutions  of  the  approximate  relation  (  ref. 3  ) 


2G1*+2aehiG-3GGi»+[G»-2ae(h1+h»)lGi  *  0 

(C-5) 

viz  Ga  =  G/2  -  p. sin((/3 ) 

(C-6) 

/a,(hi+h,)+<G/2)* 
where  p  =  2J-* - ■ - j - 

(C-7) 

.  _ i  2a  G(hs-hi) 
t  =  »iu  1  E — p - 

(C-8) 

and  where  ha  and  ha  are  the  antenna  and  target  height 
which  the  slant  range  is  calculated 

respective?'  iron 

R  -  J( h,-h,)»+4(ae+h1)(ae+h,).ain*I«!1+G,)/2a,) 

(C-9) 

Ki  (  Ra  )  is  obtained  by  substituting  ha*0  (  hi*0  )  and 
solutions  for  the  angles  6a,  6a  aud  y  are 

Ga*6  (  G»«6  )  .The 

_l  2a.(h»-ha)+h,»-h»»-R* 

"  ‘ia  - 2(.e^)R - 

(C-10) 

2a«ha+ha*+Ra* 

'  *in  2(a*+hj)Ri  ' 

(C-U) 

.  2a.ha+ha»-R*» 

*  -  »ia  2a^lT - 

(C-12) 

The  divergence  factor  (D)  accounts  for  the  faot  that  the  ray 
reflected  from  the  curved  earth's  surface  diverges  more  rapidly  than  If 
the  earth  was  flat.  For  modarata  gracing  angles  (  0<G<Go  ) 

I  a,fl. alny.coar 

D  “  /< 2GiGs/cosT+a#e.siny) <l+h»/a#) <l+h,/a#)  (C"1S) 


.../C-J 
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The  eiectroragnetic  contribution  to  the  phase  difference  (  ) 

and  reflectivity  (p)  is  obtained  from  the  coaplex  reflection  coefficient 


r  -  pe" 


(C-14) 


These  paraaeters  are  deterained  by  the  grazing  angle  (y)  and  the  coaiplex 
dielectric  constant  of  sea  water. 


(C-15) 


=  e!-j60Xo 


(C-16) 


where  a!  and  ej  are  the  real  and  iaaginary  parts  of  ec,  X  it  the 
wavelength  and  a  is  the  total  sea  water  conductivity.  Tbe  complex 
reflectivity  for  horizontal  and  vertical  polarizations  are 


siny  ~  jae- cos»y 
^h  =  siny  +  Ji 

c“COS1y” 

=  eesiny  -  ^e-cosy 
' v  e'siny  +  ^aj-cosiy 


(C-17) 


(C-18) 


The  real  and  iaaginary  parts  of  ac  are  evaluated  from  the  expressions 


(6,-4. 8)  ll+(Xt/X)  1~B(nxlO~z)  ] 

'*  “  4*8  +  1+2 (X,/X)  ^_o(nxl0”’^)  +(X,/X) 

(a  ,-4. 8)  (X,/ X) 1-a _  +  18CV, 

**  “  l+2(X./X)^-a(nxlO_2)  +  (X-/X)2^_a^  f 


(C-19) 


(  C-20) 


where  ai-0.02,  f  is  in  GHz,  X,  it  the  wavelength  corresponding  to  the 
dielectric  relaxation  frequency,  a,  it  the  static  dielectric  constant  and 
Oj  is  the  ionic  conductivity  in  aho/aetre.  Variations  in  salinity  and 
surface  temperature  are  taken  into  account  through  the  empirical 
expressions  (  reference  24,  25  ) 


a,  -  87.8  -  15. 3N  -  0.363T 


(C-21) 


100X,  -  3.38  -  0.11T  +  0.00147T*  +  0.0173TO  -  0.52N  metres  (C-22) 


oj  -  5N  +  0.12TN  +  0.04T 


(C-2S) 


in  which  T  is  the  temperat ore  in  degrees  Celsius  snd  N  is  the  electrolyte 
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normality  (  equ iva lent */ 1 it re  ).  The  relevant  temperature  it  the  water 
‘skin'  temperature,  and  the  default  value  in  the  model  ia  the  aea  level 
air  temperature.  The  default  value  for  salinity  ia  N=0.6  (3.4%). 

Equations  C-19  and  C-20  include  a  parameter  a  which  accounts  for 
the  non-Lorentzian  spectral  distribution  due  to  the  spread  in  relaxation 
frequencies  in  the  region  of  15  GHz.  Contributions  to  e*  from  water  can 
be  ignored  at  UHF  frequencies  where  it  is  dominated  by  the  18o£/fc60Xo^ 
term  as  shown  in  figure  C-2  in  which  e±  and  t%  are  plotted  up  to  100  GHz 
for  T=0,  15  and  30SC.  At  higher  temperatures  the  dielectric  relaxation 
processes  are  faster  (  l/t»20GHz  )  resulting  in  strong  temperature 
dependence  of  and  ej  in  the  frequency  range  of  interest  shown  in  the 
expanded  scale  of  figure  C-3 . 


Dielectric  subroutines,  eaq>loying  the  relations  described  above, 
were  used  to  plot  the  phase  and  magnitude  of  the  reflectivity  in  figures 
C-4  and  C-5  (  the  apparent  2n  phase  difference  at  normal  incidence  is  due 
to  the  sign  convention  for  0).  Results  for  1  and  35  GHz  display 
significant  temperature  dependence  of  p  and  0  for  vertical  polarization. 

The  dielectric  reflectivity  (p)  is  multiplied  by  the  roughness 
factor  (r)  to  obtain  the  overall  reflectivity  in  equation  C-3.  Reference 
1,  using  data  for  vertical  coverage  diagrams  for  several  sea  states, 
stresses  the  importance  of  a  realistic  model  for  the  roughness  factor. 
Although  the  literature  provides  £  prior  1  models  for  r  (  refs.  2,  26  ) 
these  are  in  general  not  in  good  agreement  with  experimental  data.  The 
mode]  uses  a  semi-empirical  relation 


I 


r 


e"2** 
e-l. 2732s 


s <0.6366 
otherwise 


(C-24) 


where  s  «  nH1^jSiny/2X  (C-2S) 

in  which  H^g  ia  the  significant  wave  height,  which  is  approximately 
related  to  the  Douglas  aea  state  Sq  by 

®l/3  *  °-5  V  <C"2<) 


for  ia  feet. 
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The  Diffraction  Ret  ion 

Mechanises  for  propagation  oyer  the  horizon  include  diffraction, 
ducting,  troposcatter  and  ionospheric  reflection.  In  the  absence  of 
•nouoloos  propagation  and  at  the  frequencies  and  ranges  of  interest  for 
radar  predictions,  diffraction  is  the  most  important.  The  model  sssur.es 
that  for  distances  greater  than  1.05  tines  the  horizon  range  the  single 
node  approximation  to  the  diffraction  field  is  adequate  (  ref. 3  ).  Using 
standard  nonenclature 


F  =  f(01)^U(Z1).D(21).V(X)' 


(C-27) 


where  Zt,  Z2  and  X  are  the  antenna  and  target  heights  and  target  range 
(  in  natural  units  )  respectively,  and  02  1*  the  elevation  angle  to  the 
horizon.  The  natural  unit  of  range  is  (  ref  27  ) 


L  =  (ae»X/it) 


1/3 


(C-28) 


and  the  natural  unit  of  height  is 
B  •*■%< aeX*/n*)1/3 

The  attenuation  factor  in  dB  is  then 
V(X)  <=  10. 99+10  log  2,X-17.55X 


(C-29) 


(C-30) 


For  the  height  gain  factors  the  aiodel  uses  the  approximate  expressions 
given  by  Blake  (  ref. 2  ) 

U(Z)  -  201ogt(Z  Z<0.6 

-  -4.3+51.04 [logs* (Z/0.6) ] 7,4  0.6fZ<l 

-  19.85(Z°-47-0.9)  otherwise  (C-31) 

The  Intermediate  Ration 

Near  the  horiion  (  y-0  )  the  calculated  divergence  factor 

approaches  zero  and  the  ray  theory  express  lone  would  predict  F~l,  whereat 
F~0  near  the  horizon.  In  the  model  ray  theory  is  only  used  up  to  the 
critical  ground  range 
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r 


G«  =  »etn-27c-*iB  1(»eeo*Yc/(ht+»e))-*lB_1(*eeo»ifc/(k»*«#))) 

(C-32) 

where  yc  is  the  critical  (rating  angle 

yc  =  tan_1(X/2nae)1/3  (C-S3) 

Between  this  range  and  1.05  tines  the  horiton  range  it  is  assumed  that  F 
decays  exponentially  (  i.e.  linearly  in  decibels  ).  This  simple 
interpolation  is  used  since  the  predicted  detection  range  in  the 
intermediate  region  is  likely  to  be  rnch  more  sensitive  to  a  realistic 
model  for  target  scintillation  than  the  interpolation  technique. 
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ANNEX  D 

ANTENNA  PATTERN  FPNCTION 

The  power  gain  (G)  describes  the  distribution  in  elevation  and 
aziamth  of  the  power  available  at  the  antenna.  In  the  Model  this  appears 
as  the  product  of  the  power  along  boresight  (G,)  and  the  normalized 
antenna  pattern  function 

G<©.0)  *  G0E* <©,#) 

“  <*«Ei (©) .E» (•)  (D-l) 


where  G« 


E»(©,#)sin©  dG  d« 


for  an  antenna  of  effiency  factor  q.  The  required  Inputs  to  the  antenna 
subroutines  are  (i)  general  type  of  antenna  (  pencil,  etc.  ).  (ii) 
horizontal  and  vertical  3-dB  beamwidths  (0H,0V),  (iii)  maximum  power  gain 
and  (iv)  the  level  of  the  first  vertical  sidelobe.  If  9®  and  ©^ 
are  known,  but  not  G(.  the  default  value  is 


C,  -  0.8  (4n/9HGV)  (D-3) 

E(0)  can  be  calculated  fro*  the  current  distribution  across  the 
aperture,  however  usually  only  a  few  far  field  parameters  are  available 
fro*  radar  manuals.  These  parameters  arc  used  in  the  model  to  synthesize 
two  general  classes  of  antenna  (i)  modified  (sin  u)/u  and  ( i 1)  cosecant 
squared.  These  patterns  form  the  basis  functions  for  other  patterns.  If 
the  effects  of  aperture  blockage  are  important,  the  pattern  may  be 
synthesized  from  a  linear  combination  at  described  in  reference  1. 


This  pattern  results  from  a  one  parameter  aperture  distribution 


of  the  form 


A(z)  -  jji  J.<iwB,/l-4z*/d*'  )  (D-4) 

where  J,  is  the  Bessel  function  of  the  first  kind  and  z  is  in  the  range 
14/1.  The  resulting  pattern  la 


.../D-a 


-  *e 


ivr- 


E(0) 


sin  itJut-B1 
n  /  u*-B*1 


(D-5) 


Equation  D— 5  is  valid  for  real  or  imaginary  arguments  and  is  equivalent  to 


sinh  n/B*-n*' 
E(®)  “  nJ  B»-n*' 


(D-6) 


B  is  calculated  from  the  sidelohe  level  by  iteratively  solving  the 
relat ion 


SL  =  (SL)«  (D-7) 

where  (SL)0  is  the  sidelobe  voltage  of  the  (sin  u)/u  pattern  (~4.603) .  To 
normalize  D~5  and  D~6  the  RHS  is  divided  by  the  noraalized  sidelobe  ratio 
(  SL/4.603).  This  ratio  also  deteraines  the  aperture  efficiency  through 
the  relation 


Jl  sinh  n7B»-(p»/2)? 
n/B*-(Pe/2)*' 


(D-8) 


with 


P,  -  (d/X)  O3®  (D-9) 

Equations  D-7  and  D-8  are  solved  iteratively  with  the  initial  estiaates 

given  by  the  polynoaial  relations 

B*  -  -0. 6604597+0. 03086396*+0. 001323133** 

+9. 498675X10~6**-1. 040798X10 _7**  (D-10) 

and  p»  -  0. 58114+0. 023468*+3 .6432X10-5** 

-6.7425X10"8**+7,9701I10“8*4  (D-ll) 

where  z  is  the  level  of  the  first  sidelobe  is  dB.  The  solutions  for  B* 
and  p(  are  plotted  in  figure  D-l.  Representative  plots  of  E(0)  using 
equations  D-5  to  D-8  are  shown  in  figures  D“2  to  D~4. 

Cosecant  Squared 


There  antennae  are  designed  to  detect,  with  constant  return. 


targets  at  constant  altitude  or  displacement  off  the  boresight  (  if  the 
tilt  is  non-zero  ).  The  corresponding  pattern  is 


where  Eo(6)  is  the  pencil  pattern  described  above.  The  gain  of  the 
cosecant  squared  antenna  is  less  than  the  value  determined  by  equation 
D-3  by  a  factor 


(D-13) 


The  default  value  of  G  for  cosecant  squared  antennae  is  3dB  less  than 
equation  D-3  (  i.e.  2©BBX>>6^db  is  assuaed.  )  Plots  of 
seaborne  and  airborne  cosecant-squared  patterns  are  shown  in  fignres  D-5 
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Laval  Of  First  Sidalobe  (  dB  ) 

Figure  D-l.  Input  Farasetsrs  For  Modified  (sin  u)/u  Pot  tarn 


Anglo  Off  Boraelght  (  dogrooo  ) 


Figure  D-2.  Modified  (oin  u)/u  Antonno  Pot  tom  Function 
Lovel  of  Firet  Sidelobo  io  10  dB  doon 
Half -pow  Boaooidth  io  1  degree 


Anglo  Off  Boreelc^it  (  dagrw  ) 


Figure  D-3L  Modified  (tin  u)/u  Antenna. Pattern  Function 
Level  of  Firet  Sidelobe  is  13. 26  dB  doen 
Hoi  f -power  Boamldth  le  1  degree 


Anglo  Off  Borooi^it  (  dogrona  ) 

Figure  0-4.  Modified  Coin  u)/u  Antamo  Pot  tom  Function 
Laval  of  Firot  Sidolabo  io  17.6  dB  doon 
Half-povar  Boaovidth  to  1  dogroo 


Antenna  Pattern  Fi 


Antenna  Pattern  Function  <  dB 
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E-l 


The  antenna  noiae  temperature  accounts  (or  the  noiae  fro*  th« 
black  body,  cosmic  and  nan  made  sources.  If  the  antenna  senses  n  noiae 
sources  of  brightness  teaiperature  then 

n  r2n  <«i 

Ta  -  ^2  Tbi(O.#)G(e,0)sin9  d©  d#  (E-l) 

isj/o  ■'o 

where  G(©,0)  is  the  antenna  gain  as  a  function  of  elevation  and  azimuth 
Tb^(9,0)  is  the  noise  temperature  as  a  function  of  elevation  and 

azinuth. 


If  the  angle  dependence  T^^C©,#)  of  each  noise  teaiperature 
source  is  known  then  equation  E-l  can  be  solved  numerically.  However,  in 
this  model  a  simplified  approach  is  adopted,  where  the  source  is 
considered  to  be  contained  in  a  solid  angle  of  average  gain  G  and  average 
brightness  temperature  and 


T.  'X^bi  <*-*> 

i 

where  k^  is  the  gain  weighted  solid  angle  divided  by  4n.  For  reasons 
discussed  later  contributions  fro*  the  sea  are  ignored,  and  for 
atmospheric  and  cosmic  sources  k^sO.9  the  following  values  are  assumed. 


k  T 

(E-3) 

CO$  COf 

lA 

ktun^"sun 

(B-4) 

kat*Tatn 

-  290  <1-1/La) 

(E-5) 

where  f  is  the  frequency  in  hertz,  the  resultant  temperatures  are  in 
Kelvin  and  LA  is  the  total  one-way  atmospheric  loss  along  the  antenna 
bores ight. 

The  total  one-way  atmospheric  loss  along  the  boresight  ia 
calculated  ia  a  manner  similar  to  that  discussed  ia  Annex  B  for  antenna 
to  target  absorption. 


E-2 


-  D-aj 


where  i*  the  total  one-way  atmospheric  loaa  due  to  constituent  j. 
However,  in  this  case  it  is  necessary  to  irtegrate  along  the  raypath  from 
the  antenna  to  infinity.  There  are  two  cases  to  be  considered  (i)  where 
the  tilt  angle  of  the  antenna  is  greater  than  the  elevation  angle  of  the 
horizon  (  and,  therefore,  the  ray  points  directly  at  the  sky  )  and  (ii) 
where  the  tilt  angle  is  less  than  the  elevation  angle  of  the  horizcn  ( 
and  the  ray  is  reflected  fron  the  sea  surface  ). 


"  Y joj  f(t,d* 

Jo 

=  Y jo  J  f(s)ds  +  I  g ( s)ds 

*  rt  •  n 


when 


otherwise 


where  f(s) 


is*-*2s(ae+h1)ain0t  +  <ae+h1)1  -ae 

T 


g(s)  =  e 


(E-10) 


T#  =  cos  1I(l+h1/ae)cos®t] 

-  grazing  angle  of  boresigbt  ray 


(E-U) 


R  = V  «e*+4ae (Sj+hi) sin* [<0t+y«)/2J 


(E-12) 


•t  is  the  antenna  tilt  angle 
I)*  is  the  antenna  height 

ae  is  the  earth's  effective  radins  (  see  Annex  A  ) 


8^  =  sin 


(E-13) 


«  the  elevation  angle  of  the  horizon 
Hj  is  a  constant  for  constituent  j 

yQj  is  the  sea  level  absorption  for  constituent  j  (  see  Annex  B  ) 


The  integrals  in  equations  E-7  and  E-8  do  not  lend  theaselves  to 
analytical  solution  and  numerical  Methods  described  at  annex  B  are  used. 
The  integrals  fro*  zero  to  infinity  are  approxiaated  by  integrals  f roa 
zero  to  an  effective  ataospheric  liait  I(,  which  is  the  range  that 
coincides  with  an  assuaed  ataospheric  thickness  h#  (  150  ka  in  the 
node 1  ) . 


(E-14) 


J 


E-3 


Re  -\( h1+*e)*-(hi+«e)*co*#7  -  (ae+hi)sin#t 


This  enables  equations  E-7  and  E-8  to  be  approximated  by 


LAj 


Tjo8e. 


,wif<*i-8e) 


•hen 


(E-15) 


=  T 


jo 


8^*if(zi.R)  + 

i=l  1=1 


otherwise  (E-16) 


where  and  w^  are  tabulated  in  column  k=0  of  table  25.8  of  reference  15 
for  various  n.  In  this  model  n*=8  is  used,  enabling  the  same  subroutine 
to  be  used  as  is  nsed  to  solve  equation  B-6. 


Equations  E-3  to  E-5  are  based  on  the  assumption  that  all  cosmic 
and  black  body  radiation  enters  the  antenna  directly  or  by  reflection 
from  a  perfectly  reflecting  sea.  With  airborne  radar  the  boresight  may 
be  pointing  into  the  sea  at  near  the  quasi-Brewster  angle.  The 
brightness  temperature  of  the  sea  is 

Tb(«.)  ■Tt(...)<1-B<®»*  (E'17) 

where  Tt(sea)  the  thermodynamic  temperature  and 
R(O)  is  the  Fresnel  reflection  coefficient. 


Using  the  subroutine  described  in  annex  C  for  H(8),  has  been 
plotted  in  figure  E-l.  At  the  maximum  for  (  5-10*  )  the  error  in  T# 
using  only  equations  E-3  to  E-5  is  in  the  order  of  a  few  degrees  and  is 
therefore  neglected. 

In  the  cates  where  Tg  does  not  contribute  signif ieantly  to  the 
total  noise  temperature,  such  as  receivers  with  NF  >>  5dB,  approximate 
formulae  may  be  used.  Reference  1  gives  approximate  expressions  for 
ground  based  radars.  Equations  E-3  to  E-5  are  used  with 

LA  m  y,  [3. 5+107(1-^) 14 1  (B-18) 

where  ®  is  the  elevation  angle  and  y0  *•  the  sea  level  attenuation  in 
dB/n.nile. 
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ANNEX  F 

CLUTTER 


Sea  Chttjf 


The  radar  equation  for  tea  clntter  say  be  written  at 


Pc  *  (4ni  >R< 


(F-l) 


where  f (6)  is  the  antenna  pattern  function  for  the  ray  to  the  clntter 
cell  which  makes  grazing  angle  r  at  slant  range  R.  For  an  irregnlar  sea 
surface  more  than  one  ray  is  possible  so  that  equation  F-l  differs  from 
the  normal  radar  equation  due  to  the  implicit  inclusion  of  multipath 
effects  in  the  clutter  cross  section  o(R,y). 


Clutter  return  depends  on  the  area  of  sea  illuminated  in  the 
minimum  time  resolution  possible.  For  a  matched-filter  receiver  with 
rectangular  spectral  envelope  the  corresponding  range  resolution  (AR)  is 

AR  =  c/2B  (F-2) 


for  receiver  bandwidth  B.  For  simple  uncoded  pulses 
AR  *  ct/2 


(F-3) 


where  c  it  the  speed  of  light  and  T  is  the  pulse  width.  For  pulse 
compression  radar  the  time -bandwidth  product  of  the  transmitted  pulse 
equals  the  pulse  compression  ratio  so  that  t  in  equation  F-3  can  also  be 
interpreted  as  the  compressed  pulse  width. 

Unlike  point  scatt erers  the  clutter  normally  ’fills'  the  beam  in 
azimuth.  It  is  customary  to  assume  an  effective  beamwidth  (  9®  *  *3dB  * 
outside  which  the  clutter  return  is  negligible.  This  introduces  an 
overestimate  of  ail. 3  dB  in  the  calculated  clutter  return  for  the  types  of 
antenna  patterns  discussed  in  Annex  0.  To  compensate  for  this  a  smaller 
effective  horizontal  beamwidth  is  used  la  the  calculations. 

e8  -  o.75  e3dB  (f-4) 


.../ F-* 


The  clutter  cross  section  is  then 


o(R.y)  =  o*  A  (F-5) 
where  A  =  0.375  R9§dBCTc0np • *ec* 

is  the  'illuminated  area'  and  c*  is  the  clutter  cross-section  per  unit 
a  rea . 

There  are  numerous  models  for  c*  in  the  literature  (  refs  2-6.28  ) 
however  these  do  not  in  general  give  good  agreement  over  the  wide  range 
of  frequencies,  grazing  angle  and  pulse  widths  required  for  this  model. 
It  was  considered  more  appropriate  to  develop  a  general  semi-empirical 
model  with  frequency  (f),  grazing  angle  (y)  and  sea  state  (S)  as  the  only 
input  parameters.  The  desired  range  of  validity  was  0.4  GHz<f<35  GHz. 
00.25m*  and  0.1®<y<10*  and  0<S<6.  Nathan  son  and  others  (  refs  6,  28  )  have 
reported  a*  over  this  range  and  these  data  were  chosen  to  fit  the 
semi-empirical  parameters. 

The  grazing  angle  dependence  of  scatter  from  the  sea  surface  can 
be  divided  into  three  regions.  At  near  normal  incidence  (y  *=  90*)  the 
scatter  is  specular  and  is  greater  for  smooth  seas  than  for  rough.  At 
intermediate  grazing  angles  the  backscatter  is  only  weakly  dependent  on 
surface  roughness  ('the  plateau  region').  Below  some  critical  grazing 
angle  (yc) ,  which  depends  on  the  frequency  and  the  surface  roughness,  the 
scstter  is  greater  for  rough  seas  than  for  smooth.  The  model  is  concerned 
with  the  region  around  yc. 

At  sea  state  3,  microwave  backscatter  Is  strongly  correlated  with 
the  appearance  of  whitecaps.  This  enables  yQ  to  be  estimated  by  analogy 
with  backscatter  from  a  target  at  a  constant  height  and  consideration  of 
multipath  effects  (  discussed  at  Annex  C  ).  Plots  of  the  pattern 
propagation  factor  (F)  for  an  array  of  scattarars  at  a  haight  of  1.5 
metres  above  sea  level  display  typical  interference  lobing,  with  the 
loaest  angle  lobe  at  a  grating  angle  near  the  limit  of  ray  theory 
applicability  (Annex  C) .  The  lowest  lobe  maxima  for  1  GHz  and  10  (Six  are 
observed  at  yc  -  3*  and  0,3*  respectively,  A  linear  interpolation  is 
assumed  for  ail  frequencies  at  sea  state  3 


Tc  * 


iMtj  *•**••• 


<F-« 


J 
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The  tea  state  dependence  of  yc  is  complex,  since  the  relative 
importance  of  physical  scatterers  (  swell,  sea,  capillary  waves  and  foam) 
changes  with  both  freqoency  and  sea  state.  In  the  absence  of  detailed 
experimental  data  near  yc  it  is  assumed  that  yc  is  inversely  proportional 
to  sea  state.  Fitting  Nathanson's  horizontal  polarisation  data  for  t>0.5 
ps,  0.1<y<10#,  and  an  average  of  npwind/ downwind  data  and  sea  state  3  : 

oMdB/mVm1)  =  -159  +  31.61og(f)  +  19.91og(Yc) 

+  1701og(r/Yc)/fM 
=  -159  +  19.91og(T)  +  31.61og(f) 
for  f  in  MHz  and  y  in  degrees. 

Data  for  vertical  polarization  does  not  exhibit  a  low  angle 
plateau  region  and  o'  (dB/n'/n* )  vs  log  y  is  fairly  linear  over  the  region  of 
interest.  For  sea  state  3  and  vertical  polarization,  the  fit  of 
Nathanson's  data  yields  : 

o'  (dB/m'/m*)  =  -30  +  410f-*'log(y/10)  +  J.51og(f/50000)  log(y/y]f) 


for  y>y^ 

-  -30  + 

410f-'*log  (y/10) 

otherwise 

(F-8) 

The  critical  angle  for 

vertical  polariration 

<Yc>  appears 

to  be 

insensitive  to  wavelength  and  wave  height  and  an  arbitary  valne  of  1*  has 
been  chosen. 

Variations  in  the  wave  height  affect  not  only  the  critical  angle 
yc,  bnt  also  flnctnations  in  numbers  and  positions  of  scatterers.  This 
contribution  appears  to  vary  as  the  fonrth  power  of  sea  state  for 
horizontal  polarization 

o'(S)  -o'(S-3)  +  401og<S/3>  (F-9) 

This  simple  behsvionr  is  not  seen  for  vertical  polarisation,  where  »• 
appears  to  vary  by  abont  4.3  dB  per  sea  state  below  sea  state  3.  Above 
sea  state  3,  where  clntter  is  strongly  correlated  with  the  appearance  of 
white  caps,  and  shadowing  oecvrs  for  low  y,  o'  seems  to  vary  by  only  3.5 
dB  per  sea  state. 

Fignres  F-l  to  F-6  show  computed  o*  for  typical  radar  and 


for  y<yc 

otherwise  (F-7) 


. . . /F-4 


F-4 
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environmental  parameter*.  The  overall  level  of  agreement  is  +  5  dB  which 
is  comparable  with  the  accuracy  of  the  data.  Wind  aspect  dependence  of  c° 
is  not  included  in  the  model  since  the  only  data,  which  is  for  larger 
grazing  angles  (10  -  60  degrees),  does  not  display  a  well  behaved  relation 
betweci  o°  and  aspect  angle  (  ref  28  ).  Most  models  assume  some  circular 
function,  such  as  a  cosine,  to  predict  the  wind  aspect  dependence  of  c°. 
From  the  available  data  it  ia  not  clear  that  this  is  valid  for  the  wide 
range  of  giazing  angles  and  frequencies  of  interest.  This  model  therefore 
accounts  for  wind  aspect  only  as  an  additional  nncertainty  of  .+5  dB  in  o° 
for  simulation  purposes. 

Non-integer  sea  states  are  allowed  in  the  model  it  order  to 
facilitate  small  variations  in  wave  height  for  simulation  purposes.  Data 
for  nominal  sea  state  0  in  reference  6  was  used  to  fit  S<=2/3  in  the  plots 
for  horizontal  polarization,  consistent  with  Kathanson’s  comments  about 
minimum  wave  height. 

Amplitude  and  Temporal  Statistics 

The  previous  sub-section  describes  the  calculation  of  median 
clutter  return  for  UHF  to  microwave  frequencies  and  low  grazing  angles. 
To  calculate  probability  of  detection  from  signal-  io-clutter  ratios  a 
model  is  also  required  for  the  amplitude  and  temporal  statistics  of  the 
clutter. 


This  model  assume*  that  the  clutter  resolution  cell  is 
sufficiently  large  for  a  large  number  of  independent  scattarers  of  random 
phase  to  be  detected.  Under  these  conditions  the  clutter  ia  Rayleigh 
distributed.  A  minimum  area  of  104  m*  seams  reasonable  at  I-band  and 
corresponds  to  r>0.2  (is  (  B<5  MHz  )  and  ®®>1*  for  typical  ranges.  This 
range  is  confirmed  by  Nathanson's  short-pulse  I-bamd  distribution 
(  ref  6  )  which  shows  that  the  ratio  of  «•  to  its  standard  deviation 
is  near  unity  for  horizontal  polarization  only  for  t>0.1|is.  Vertical 
clutter  seem*  less  spikey  even  at  vary  small  puls*  widths,  and  tb* 
Rayleigh  distribution  is  probably  adequate  for  all  puls*  widths  of 
interest. 

I  A  model  of  temporal  statistic*  it  required  for  estimation  of  some 

I  clutter  reduction  techniques.  Within  tha  range  of  validity  of  the  Rayleigh 


F-5 

amplitude  distribution,  the  clatter  appears  tc  decorrelate  on  a 
Billisecond  timescale.  This  is  not  trne  for  very  short  pulses  where  the 
deceit e lat ion  tiae  say  be  of  the  order  of  several  seconds.  Under  these 
conditions  clutter  returns  at  different  frequencies  may  not  be 
deconelated,  so  that  the  criterion  Af  >>  1/t  for  frequency  agility 

decorrelation  may  not  be  valid.  This  subject  requires  deeper  investigation 
before  short-pnlse  decorrelation  techniques  can  be  modelled  adequately. 

In  summary  the  sea  clutter  model  used  in  the  programs  is  expected 
to  produce  realistic  probabilities  of  detection  for  r>0.2  ps,  0.4<f<35GHz. 
sea  states  0  tc  6,  0.1°<y<10#,  horizontal  or  vertical  polarization  for 
standard  pulse  radars. 

Volume  Clutter 

The  volume  clutter  cross  section  (o)  is 

o  -  oyV  «  cv  •^R*6®9VAR  (F-10) 

where  R  is  the  slant  range,  9®  and  6^  are  the  effective  beamwidtbs  (~w/4 
times  the  3dE  beamwidth)  and  the  range  resolution  AR  is 

AR  =  c/2B  (F-ll) 

for  matched  filter  of  bandwidth  B. 


The  specific  clutter  cross  section  (ey)  is  strongly  frequency 
dependent.  At  UHF  and  low  microwave  frequencies  (<10GBz)  Rayleigh 
scattering  behaviour  applies  and,  ignoring  multiple  scattering  effects 


*V  -V  -frltl'Z 


(F-12) 


where  sc  is  the  complex  dielectric  constant  of  water,  dj  is  the  diameter 
of  the  ith  scatterer  and  the  summation  if  over  all  particles  in  the  clutter 
cell.  Being  the  formulae  for  se  (  for  pure  water  >  described  in  annex  C. 


Ill*  -  0.93  ±  0.01 


(F-13) 
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over  tbe  range  of  frequencies  and  temperatures  of  interest.  An  empirical 
expression  for  the  reflectivity  (Z)  is  (  ref  4  ) 


Z  =  2X10-16  I1'6 


(P-14) 


where  r  is  the  rainfall  rate  in  mn/honr.  Combining  equations  F-12  to  F-14 


(F-15) 


Snow 


Snow  can  be  treated  in  a  similar  manner  dne  to  a  fortuitous 
cancellation  of  small  |K I *  and  large  Z-r  coefficient: 

III*  a  C.20  (F-l<) 

Z  a  lO-15  r1-6  (F~17) 


where  r  is  interpreted  as  the  snowfall  rate  after  the  snow  has  aelted 
(  ref  4  ) . 
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Flgura  F-4.  Soo  Cluttar  Grow  Motion  at  900  MHz 
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ANNEX  G 

PROBABILITY  OF  DETECTION 

Probability  of  detection  is  calculated  in  tbe  model  for  fixed  and 
adaptive  threshold  (CFAR)  detectors.  Rayleigh  amplitude  statistics  are 
assumed  for  the  reciever  noise  and  clutter,  while  the  target  can  be 
non-fluctuating  or  of  a  generalized  Chi-square  distribution. 

The  general  method  for  calculating  paint  probability  has  been 
described  by  Karcum  (  ref  Id  ),  Swerling  (  ref  17  )  and  Rubin  and  Di 
Franco  (  ref  18  ) .  Recursive  algorithms  are  similar  to  those  described 
elswhere  in  the  literature  (  refs  19,  20  )  with  additional  scaling 
algorithms  to  reduce  the  dependence  of  the  results  on  the  numerical 
precision  and/or  dynamic  range  of  the  computer  in  use. 

Kon-coherent  Detection 

Assuming  Rayleigh  noise/clutter  statistics,  the  distribution  of 
the  integrated  signal-plus-noise-plus-clutter  variate  (v)  for  a  given 
integrated  signal-to-noise-plus-clutter  ratio  (a),  where  v  is  normalised 
with  respect  to  the  average  noise  level  is  given  by 


fN<vlx)  “  f**v,H,'1»',’+‘) 


where  N  *  the  number  of  pulses  incoherently  integrated 
a  *  Nx 

x  »  the  single  pulse  signal-to-noise-plus-eletter  ratio 


The  probability  that  v  will  exceed  a  fixed  threshold  (T)  la  given 


P(v>Tla) 


*  j" 

•  £■*•  rv 


For  a  non-flaataatlng  target  It  la  sufficient  to  replace  a  with  the  mean 


.../•-a 


1 


integrated  signal-to-noise-plus-clutter  ratio  (Z)  and  evaluate  this 
equat ion. 

To  determine  this  probability  for  a  fluctuating  target  it  it 
necessary  to  integrate  over  the  target  distribution  w(z,  Z)  and 


=  I  i w(z,Z) 

*  n 


P(v>Ylz)  dz 


In  this  paper  w(z,  Z)  is  assumed  to  be  in  the  form  of  the  Chi-square 
distribution  and 


'•Z)  c  HR  *~frfT 


where  K  is  a  distribution  parameter  selected  as  required  to  synthesize 
actual  target  fluctuations.  K  may  only  take  positive  values  and  is  of  the 
order  of  the  number  of  independently  fading  scatterers.  The  four  Swerling 
cases  may  be  evaluated  by  substituting  K  *  1,  N,  2,  and  2N  respectively. 
Substituting  equation  G-4  into  equation  G-3  and  performing  the  necessary 
integration  gives 


The  probability  of  false  alarm  (Pfl)  la  the  probability  that  v  will 
exceed  T  in  the  absence  of  signal  and  may  be  found  by  substituting  Z  -  0 


“fa  •  rT* 


Constant  False  Alarm  Kate  ICFAJU 

In  the  case  of  CFAR  detection  tha  fixed  threshold  T  is  replaeed 
by  a  variable  threshold  (y)  which  is  derived  from  the  sum  (n)  of  R 
independent  Rayleigh  noise-plue-clutter  samples  of  unit  average  power. 
Therefore  u  has  tha  distribution 

‘  fl'-l?  !*  <e"7) 


Substituting  y  “  u/T,  whare  T  is  a  calibrating  factor,  and  performing 
the  necessary  integration  gives  the  overall  probability  of  v  exceeding  T 


.../•-I 


Jt 


G-3 
N+a-l 


n  -  e  jflia)  n  mrnitiS  f*HT 


(G-g) 


As  before  the  probability  of  false  slant  aay  be  found  by  substituting 
Z  =  0  and 


N-l 


-u  -E'atFiHfiW'W 


(G-9) 


Alaorithrs  for  P . 


Let 

F(n)  *=  for  non-fluctuating  targets 

“  nifri]  (i$zV  $z)  for  fluctuating  targets  (G-10) 

vst  -Y 

G(a)  «  for  fixed  threshold  detectors 

vK/ 


■Itl-iH  (t+rf  ,or  t**®  d®t«etor* 


(G-ll) 


so  that 


N-l 

Pj,  •  JEG(m)  for  cases 

■*o 

(G-12) 

Now, 

]EZp(®)  *  5Z  g(o)  “i 

n*>o  n>=o 

(G-13) 

and 

-  nan 

P4  -  H  F(n)  ZZ  G(«) 

(G-14) 

n«o 


Expanding  this  and  regrouping  the  soamatious  gives 

N-l  •  n-N 

tA  -  ZIG(a)  +  216(b)  [1  -  CF(n)] 
n-o  n»N  n»o 


£Xfl(«)  +  £lG(u)Cl  -  Z?F(a)J  ♦  C  G(b>[1  -3E?F(n>] 

■*o  sq|  n“o  rt+1  a*© 


vhere  B>N 


(8-15) 


Therefore 


F4  -  |j£e<»)  ♦  |j^8(B>ll  - 


♦  a 


(8-10 


where  the  argot  ten 


.../8-4 


G-4 


®  b-N 

e  =  T-  G(a)  [1  -  ZI  F(n)  ] 
m=B+l  n=o 


B  B-N 

<  (1  -  Y  G(b))  (1  -tZ  F<n)) 
m=o  n=o 


where  B>N 


(G-17) 

(G-18) 


This  error  tern  may  be  used  to  truncate  the  infinite  summation 
when  a  specified  error  limit  is  reached.  Alternatively,  the  expressions 


N-l 

H  =  52  G(m) 
m=o 


B  B"“N 

+  T.  G(m)  [1  -  SZF(n)] 
m=N  n=o 


(G-19) 


and 


N-l 

i  =  ii;g(«i) 

m=o 


B  m-N  B  B-N 

+  ZlG(m)[l  -  2ZF(n)}  +  (1  -  5IG(b))(1  -  ZIF(n)) 
m=N  n=o  m=o  n=o 


(G-20) 


both  converge  on 

P.  »  ~  £TF(n)J  as  B->«  (fi-21) 

“  m=o  m=N  n«=o 

The  value  of  I  converges  more  rapidly  on  Pj  than  does  that  of  H  and  is 
therefore  used  in  the  model.  Summation  is  terminated  when  successive 
estimates  of  differ  by  less  than  a  specified  amount  and  the  'error 

term'  is  less  than  0.1/N.  (  This  latter  condition  is  necessary  to  prevent 
preBiature  truncation  of  the  series  for  high  signal-to-noise-plus-clutter 
ratios  ). 


Recursive  Relations 

The  terms  in  the  infinite  series  are  evaluated  recursively  using 
the  following  relations  : 

a.  F(n+1)  -  F(n) 

with  F(0)  -  e~^  for  non-fluctuating  targets  (6-22) 

b.  F(n+1)  -  ¥&  i5r  •F(n) 

with  F(0)  -  for  fluctuating  targets  (6-22) 

c.  6(m+l)  -  1^-  T$T  .6(a) 


.../6-S 


G— 5 


with  G(0)  *  (f+r)  for  CFAR  detection  (G-24) 

d.  G(m+1)  =  ^  ,G(m) 

with  G(0)  =  e-^  for  fixed  threshold  detection  (G-25) 

Evaluation  of  Threshold 

In  most  applicstions  of  the  model,  the  nser  wishes  to  calculate 
the  performance  of  the  radar  for  a  given  false  alarm  rate  or  probability 
of  false  alarm  and  it  it  therefore  necessary  to  calculate  a  valne  for  the 
threshold  as  a  function  of  the  probability  of  false  alarm. 

When  the  number  of  pulses  integrated  (N)  is  1.  the  solutions  are 

given  by 


Y  =  -  loge  (  P£t  )  for  fixed  threshold  detection  (G-26) 

and 

T  -  2  /  (  1  -  Z  )  for  CFAR  detection  (G-27) 

where  Z  =  exp  [  /  *  J 

When  N  >  1  ,  there  is  no  analytical  solution  and  the  value  of  Y 
is  determined  by  Newton’s  method.  The  relationship  between  Y  and  Pj# 
described  at  equation  G-6  does  not  satisfy  all  of  the  pxerequisites  for 
Newton's  method  and  therefore  the  following  well-behaved,  normalised 
function  is  used 

g(X)  -  loge[  f(X)  /  Pf,  ]  (0-28) 

and  succesive  iterations  for  the  solution  of  g(X)  -  0  are 

Vi  -  Tk  -rr^y,,<Tk) 

For  the  fixed  threshold  detector  T  replaces  X  and 

t-1  -Y.n 
f(Y)  “  )  \'r 
■» o 


(0-29) 

(0-90) 
. . ,/o-d 


-YyN-1 

f'<Y>  =  -  TJFTTT 


(G-3 1) 


which  yield* 


N-l  N-l 

5ZG(m)  jr  G(a) 


?K+1  “  YK  +  (Hti-l)  >lo8 


(G-3 2) 


For  CFAR  detectors  T  replaces  X  and 


f(T)  =  ^mHR-ITT  (t+T^ 


(G-33) 


'«)  -  g1a^fT^H(*)R(rhr 

D  N-l  N 

4E6(«)  -  Z>G(») 

B=0  11=1 


(G-34) 


Nnaiber  of  Pulses  Integrated:  2-D  Radar 


The  number  of  pulses  incoherently  intgerated  (N)  for  a  scanning 
radar  is  assumed  to  be 


N  „ 


(G-35) 


where  is  the  horizontal  beamwidtb, 

prf  is  the  pulse  repitition  frequency  in  Hz,  and 
4  is  the  target/bore* ight  azimuth  rate 


4  **  6rp* 


180  V.  sinO. 


(6-30 


in  which  K  is  the  target  range,  V  ia  the  target  velocity  in  range  units 
per  sec  and  6r  is  the  angle  between  the  radar-target  aais  and  the  target 
velocity  vector.  The  effect  of  the  target  ten  is  neglected  for  routine 
calculations,  since  it  is  tlgulfietat  only  for  very  short  range  crossing 
targets  with  snail  N.  In  the  preceding  sections  K  is  rounded  down  to  an 
integer  value. 

ULli&U 


Simultaneous  aaiamth  and  elevation  scanning  is  normally  performed 


7 


J 


G-7 


by  phased  array  antennae.  The  number  of  hits  per  scan,  ignoring  target 
aiotion,  is 


N 


eVprf 

TSTT^rpai 


(G-37) 


where  uy  is  the  vertical  scan  rate,  ty  is  the  vertical  scan  period. 
Typical  values  for  3-D  target  indication  radar  are  (eH.6V)*4  deg* ,  lOrpm, 
prf  -  500Hz,  =  10'deg/sec,  ty  =■  30msec  yielding  N  »  1. 


Elevation  scanning  is  often  performed  by  frequency  jumping  in 
discrete  increments  to  provide  one  pulse  per  beam  which  are  typically 
about  t  0^  apart  to  ensure  continuous  vertical  coverage.  This  means  that 
there  is  a  significant  contribution  up  to  the  2nd  or  3rd  beam  off  the 
line  of  sight  and  increases  the  effective  number  of  pulses  per  scan  by  a 
factor  of  2  or  3.  If  these  pulses  are  integrated  on  a  PPI  an  additional 
collapsing  loss  (Lc)  must  be  included. 

Lc  -  A0/0V  (G-38) 


where  A8  is  the  total  coverage. 


Coherent  Detection 


The  radar  equation  for  coherent  processing  can  be  generalised  in 
terms  of  signal-to-noise-plus-c lutter  energy  ratios 


„  (NcPtT)Ge*X»»F4 
%  m  £t,u»!)»r4l 


(G-39) 


where  t  it  the  radiated  pulse  width.  The  equivalent  numbei  of  pulses  Nc 
is  considered  for  three  detectors  i)  simple  coherent  processing,  (ii) 
pulse  doppler  and  (iii)  pulse  compressicn. 


With  simple  coherent  detection  Ne  is  the  number  of  pulaee 
integrated.  Assuming  no  post  detection  integration  (Nj»i)  the  threshold 
level  ie  given  ae 


T  -  -lo«0<Pf,>  (G-40) 

Figures  6-12  end  6-13  show  comparisons  between  pure  coherent  end 


J 
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incoherent  integretors  with  a  non-fluctuating  target.  For  a  constant  N> 
the  incoherent  processing  loss  is  only  weakly  dependent  on  the 
signal-to-noise-plns-clntter  ratio.  Also  near  nnity  signal-to-noise-plns- 
clntter  ratio  the  'integration  loss’  varies  only  by  a  few  dB  between  1 
*nd  100  pulses  and  can  be  represented  by  the  simple  empirical  expressions 
described  in  the  literature  (  ref  S  ). 

Pulse  doppler  radar  is  also  described  by  equation  G-39  with 

VNc  *  Pt*P**/Bn  =  P.v/Bn  (G"41) 

where  Bn  is  the  doppler  filter  bandwidth,  assuming  that  the  range  gate 
is  of  the  order  of  one  pulse  width. 

The  analysis  of  pulse  compression  radar  it  essentially  the  same 
as  for  a  radar  which  radiates  at  the  compressed  pulse  width  rc  with 
increased  power  PtNc  where  Nc  is  the  time-bandwidth  product  (  or  pulse 
compression  ratio  ) 

VcKc  *  Pt*  *  ptNc/Bn  <G~42> 

The  effect  of  additional  non-coherent  integration  has  already  beau 
described.  Frequency  weighting  is  usually  employed  to  reduce  the 
time/range  sidelobes  below  the  13  dB  for  unweighted  linear  FN.  The  family 
of  Taylor  weightings  is  commonly  used,  and  are  specified  by  a  single 
parameter  n .  As  n  becomes  large  the  sidelobe  reduction,  for  a  given  pulse 
broadening,  approaches  that  of  the  Dolph-Chebyshev  weighting  at  the 
expense  of  total  sidelobe  energy.  Typioal  performance  figures  are  50% 
pulse  broadening  for  25  to  45  dB  sidelobe  level  resulting  in  a 
processing  loss  of  0.5  to  1.5  dE.  The  default  value  in  the  model  is  1.5 
dB. 


Acquisition  Probabilities 

The  probability  of  detection  is  the  probability  of 
achieving  m  paints  out  of  n  successive  trials  (scans)  for  the  first  time 
on  the  Nth  trial.  Cumulative  probability  of  detection  is  the 
probability  that  the  target  haa  been  initially  detected  by  the  Nth  trial, 
and  the  complement  G^(*)  is  the  probability  that  no  initial  detection  has 


occnred  by  the  Nth  trial. 


The  previous  section  describes  the  calculation  of  probability  of 
a  paint  pj  for  a  given  scan  i.  Recursive  algorithms  for  computing  Pj^(J) , 
P£|(*)  and  *re  described  belov 

P£(JJ)  «  0  for  N  <  m 

=  1-  Qjjj(")  otherwise  (G-43) 

pn<;>  -  pn<;j  -  ps-i<;> 

-  -  «&<;>  (G-44> 

The  expansion  of  these  expressions  becomes  more  unweildy  as  more  complex 
detection  criteria  are  evaluated  and  no  generalised  m/n  probability 
algorithms  in  terms  of  pi  are  possible.  In  the  current  model  three 
speoifio  detection  criteria  are  addressed,  i)  one-out-of-one,  (ii) 
two-out-of-two  and  (iii)  two-out-of-three,  and  algorithms  axe 

1/1  Detection  Criterion 

«S<J>  -  (1  -  Pl) 

-  (1  -  PN>«N-i<J)  for  N>1  (G-45) 

2/2  Detection  Criterion 

<*!<:>  *  i 

«§<J)  -  (1  -  P»?n-1)QN-1(;)  for  N>1  <0-4«) 

2/3  Detection  Criterion 

0|<J>  -  1 
«§<•>  -i-»ip2 

«$<;>  -  (i  -  rjrjMi  -  p,>  ♦  pjU-»a>u  -  px> 

flft(j)  -  (1  -  »M>«&-1<:>  ♦  »n<*  -  *im»«  -  P*-a>«&-s<:> 

for  N>J  (0-47) 

The  development  of  similar  roenrsive  expressions  for  other  m/n  oriteria 
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using  Markov  state  diagrams  has  been  described  by  Caste!  U  (  ref  21  )  and 
Post ema  (  ref  22  ). 

Operator  Efficiency 

The  probabilities  p^  are  probabilities  of  paint  and  do  not  take 
into  account  degradation  in  operator  performance  dne  to  fatigue,  etc. 
This  results  in  an  unrealistically  narrow  distribution  in  the  detection 
range.  In  the  model  pj  is  scaled  by  an  operator  efficiency  factor  pCp 
which  represents  the  probability  that  the  operator  will  see  a  paint  given 
that  a  paint  occurs, 

Ac c,u i s  it  ion  Range 

Two  quantities  of  interest  are  i)  what  is  the  probability  that  a 
detection  occurs  on  a  given  target  run  and  (ii)  given  that  detection 
occured,  what  is  the  expected  range?  The  former  i*  s imply  the  cumulative 
probability  of  detection  . 

The  expectation  values  of  the  detection  delays  N  and  N*  are  also 
calculated  in  the  programs 


N  -  INP£ 

N*  c  5.  N*Pj!j  (G-48) 

where  P£  =  PN/JJPN 

from  which  the  aiean  and  standard  deviation  of  the  acquisition  time  can  be 
calculated 

T  -  N  fit  (G-49) 

ox  »  fit/iT*  -  N*  (G-JO) 

where  fit  is  the  time  per  scan.  Also  calculated  axe  the  aedien  delay 
(?£•(") "0.5)  and  quart ile  (P^<") "0.23,0.79)  delays.  These  are  especially 
relevant  with  multi-aiodal  distributions  produced  in  the  nnltipath 
environment . 


Probability  of  Detection 
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Figure  G-12.  Probability  of  Detection  (  Blip/Scan  Ratio  > 
Non-fluctuating  Target,  N-10 
Comparison  of  Coherent  and  Non -coherent  Integration. 
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ANNEX  B 


1I0VXKC  TABGfT  INDICATOR  (ATI) 


Returns  from  zero  doppler,  pnlse-to-pulse  correlated  targets,  are 
reduced  in  DTI  systems  by  a  series  of  n  delay  line  cancellers.  The  model 
assures  that  the  HTJ  signal  retnrn  for  an  n-delay  canceller  is  modified 
by  a  factor 


(H-l) 


where  v  is  the  range  rate  and  >  is  the  wavelength.  Optimum  target 
velocities  occur  when 


pff  "  $<*■-!>  <H"2> 

for  integer  values  of  m.  and  blind  speeds  occur  when 

p-H-'f  <®-3> 

The  signal  gain  when  averaged  over  all  target  velocities  Is 
obtained  by  integrating  equation  H-l 

5  ‘  W*r  ‘  <B~4> 

Then  detection  is  receiver  noise  limited  the  increase  in  noise  power  is 

*  W-  <H-5> 

so  that  the  signal-to-nolse  ratio  remains  constant  in  the  absence  of 
clutter,  when  averaged  over  all  target  velocities.  The  effective  number 
of  pulses  for  post-detection  integration  is  (  ref  2  ) 


eff 

for  n*l 

(H-«) 

'eff 

for  n«2 

(1-7) 

'eff 

for  »*1 

<■-$> 

where  N0  =■  N-n  end  K  i*  the  number  of  poises  input  to  the  cenceller. 


Performance  in  clutter  it  determined  by  the  pulse  repetition 
frequency  and  the  spectrum  of  the  clutter  return  W(f).  If  the  gain 
profile  i» 

G(f)  -  22n  si n2n(f|r)  <■"»> 

and 

-f*/2«  * 

W(f)  -  W„e  c  (H-10) 


the  clutter  power  at  the  output  of  the  n-stage  MTI  is 
,+prf/2 

Poc  -  W(f)G(f)  df  (H-ll) 

JT>rf/2 

~  [  W(f)G(f)  df  (H-12) 

J  — <C 

-  W.(is/n[(1.3.J....(2n-l))^n]ffc2“+1  (H-13) 

while  the  corresponding  clutter  power  entering  the  canceller  is 


pic  “  *e«0 


(H-14) 


A  measure  of  the  ability  to  tee  an  'average'  noving  target  in 
clutter  is  the  gain-weighted  cancellation  ratio  or  improvement  factor  (I) 


I  - 


(■-15) 


Most  MTI  systens  use  a  linear  limiting  aaplifier  preceding  the 
canceller,  which  ensures  that  the  residual  clutter  resembles  receiver 
noise.  An  undesirable  effect  of  limiting  is  to  reduce  the  effective 
value  n  to  n'  where  n 'wn  for  no  limiting  and  n'ar0.9  for  hard  limiting. 
Whilst  this  reduces  I  for  single-delay  MIT  by  only  a  few  dB  it  ean 
degrade  the  performance  of  multiple-delay  KIT  systems  by  tans  of 
decibels.  This  uncertainty  juatifias  the  use  of  simple  velocity  spectra 
for  tea  and  rain  o letter,  via 

•T  -  / a  for  sea  e latter 

(or)*  -  1  +  (O.OXeV  m/a  for  rain  e letter 


where  S  is  the  sea  state  ( 

5  n/sec/kc  of  altitude  ),  the 

by  equation  (  H-10  )  and 

Douglas  ),  I  is  the  wind  shear  (  up  to 

doppler  spectre  cf  surface  clutter  is  given 

®c  c  2«VA 

(B-16) 

The  variability  in  I  for  sinulation  purposes  is  accommodated  by  varying  oy 

and  n'  (  between  0.9  and  n  ) 

In  the  evaluation  of  S/N  =  Pr/ (PC+PB)  for  the  probability 
algorithms,  the  levels  at  the  output  of  the  MTI  are  asswed  to  be 

fPr)out  -  pr-e<f> 

(H-17) 

(Pnlout  Pn*® 

(H-18) 

(PcW  ^  Pc-G'1 

(H-19) 

If  only  'average'  MTI  porfornance  is  considered  G(f)  »  G  and  so 
S/N  in  the  absence  of  clutter  is  independent  of  the  MTI  gain.  For 
s inula t ion  purposes,  where  target  velocity  is  known,  it  is  important  to 
determine  whether  the  target  is  located  near  a  blind  doppler  region 
(  Equation  H-3  ).  Often  the  HTI  pulse  repetition  frequency  is  varied, on  a 
scan-to-scan  basis,  between  two  frequencies  f»  and  f*.  If  f»/f*«l+l/2» 
(  where  b  is  a  Kultiplier  of  the  first  blind  speed  )  then  for  typical  target 
speed,  average  performance  will  be  obtained  overall  at  this  target 
doppler  although  the  target  return  nay  be  severely  degrsded  on  alternate 
scans.  In  the  worst  case  a  target  nay  renain  in  a  blind  tone  even  with 
alternating  pulse  repetition  frequencies. 
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EXAMPLES  OF  PROGRAM  OPTPUT 

Figures  1-1  to  1-6  illustrate  the  graphical  output  of  the 
programs  for  a  variety  of  target  and  detection  parameters,  as  described 
in  annex  G.  All  samples  are  for  the  sane  radar  and  environmental 
parameters  which  are  summarized  in  table  1-1. 

In  figure  1-1  the  aingle  pulse  signal-to-noise-plus-clutter  ratio 
is  plotted  against  target  ground  range  for  the  interference  and 
intermediate  regions.  Under  the  non-ducting  conditions  assumed  here,  ray 
theory  remains  valid  from  zero  range  to  the  peak  of  the  lowest 
interference  maximum  (at  approx.  20  n. miles) .  The  barely  perceptible 
change  in  slope  near  20  n. miles  justifies  the  use  of  the  exponential 
interpolation  for  the  intermediate  region.  For  ranges  less  than  the 
clutter  horizon  (12.1  n. miles),  the  S/N  envelope  begins  to  be  dominated 
by  the  clutter  return  at  sea  state  3. 

The  probability  of  detection  computational  procedures  of  annex  G 
are  illustrated  in  the  remaining  plots.  In  figure  1-2  the  S/N  data  has 
been  used  to  compute  the  paint  probabilities  for  a  scan-to-scan 
fluctuating  target.  The  30%  paint  probability,  or  unity  signal-to-noise 
ratios,  are  often  used  to  produce  nominal  detection  ranges  -  from  these 
data  both  criteria  would  lead  to  a  computed  detection  range  of 
approximately  23-24  n.m.  Even  these  crude  detection  criteria  illustrate 
dramatically  the  importance  of  including  multipath  effects  since  the 
notional  detection  range  is  only  about  75%  of  the  radar  horizon  at  S  GHz. 
This  degradation  against  low  altitude  targets  is  even  greater  for  UHF 
radars  (  ref  1  ) . 

Performance  against  very  slowly  fluctuating  targets  is 
illustrated  for  the  Veinstook  target  (K  -  0.4)  in  figure  1-3.  Between  20 
and  40%  paint  probabilities  there  is  negligible  difference  from  the 
scan-to-scan  scintillation  rate  previously  discussed.  The  differenee  is, 
howevsr,  fairly  large  in  the  region  of  greatest  interest  (50%-99%). 

A  more  significant  detection  range  estimate  is  obtained  from  the 
cumulative  probability  of  initial  detection.  This  is  illustrated  in 


V, 


.../i-a 
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figure  1-4  for  the  tingle  blip  detection  criterion  end  •  range  of 
operator  efficiencies.  These  plots  illustrate  graphically  the  essential 
difference  between  applying  an  operator  loss  (typically  3dB)  and  scaling 
the  probability  by  an  operator  effiency  or  alertness  factor.  For  modest 
signs 1-to-noise  ratios  the  foraer  wonld  normally  lead  to  a  predicted 
detection  range  near  the  peak  of  the  first  aultipath  mssiaun.  By  contrast 
the  use  of  a  realistic  alertness  factor  for  a  fatigued  operator  (20%-50%) 
leads  to  the  srach  widtr  distribution  shown  in  figure  1-4  for  a  typical 
aircraft  closing  speed.  If  the  detection  criterion  is  changed  to  two 
consecutive  blips  the  degradation  is  even  aore  draaatlc.  As  shown  in 
figures  1-5  and  1-6  for  aoderately  and  slowly  fluctuating  targets  the 
nedian  detection  range  could  easily  be  less  than  50%  of  the  range 
calculated  using  simplistic  detection  range  methods. 


Table  1-1. 


Radar  and  1 


Transmitted  Frequency 
Peak  Power 

Transmitted  Pulse  Vidth 
IF  Bandwidth 

Pulse  Repetition  Frequency 

Receiver  Noise  Figure 

Plumbing  and  Miscellaneous  Losses 

Horizontal  Beaawidth 

Vertical  Beaawidth 

Level  of  First  Vertieal  SJdelobe 

Antenna  Tilt  Angle 

Antenna  Rotation  Rate 

Antenna  Height 

Polarization 

Sea  Surface  Teaperature 
Salinity 

Sea  State  (Douglae) 

Ait  Teaperature  at  Sea  Level 
Baroaetrie  Pressure  at  Sea  Level 
Huaidity 


5000  MHz 
200  kV 

1  |IS 

1  MHz 
500  Hz 
5  dB 
10  dB 
1  degree 
20  degrees 
17.6  dB 

0  degrees 
20  rpa 
100  feet 
Horizontal 

IS  degree  C 
3.4  % 

3 

13  degree  C 
1013.25  abar 
SO  % 


)  S  10  IS  20 

Tcrgnt  Ground  Ronga  ( rva.  > 

Figure  1-2.  Probcfcllity  of  Dotation  (  Bllp/Soan  Ratio  ) 
Targoti-  Soarling  Com  I 

1  o^&  Radar  Crow  Soction 
200  foot  Altitude 
Thrwhold  lo  Flxad 

Probability  of  Fdlw  Alov  lo  a 000001 


Probability  (  X  ) 


Target  Ground  Range  ( rve.  ) 


Figure  1-4.  Cumulative  Probability  of  Initial  Detection 
Detection  Criteria  ie  a  Single  Blip 
Target*-  Sear  ling  Con  I 

1  eq.e.  Rada*  Croee  Section 
200  Feet  Altitude 
480  knote  Cloeing  Rate 
Threahold  le  Fixed 

Probability  of  Falee  Alarm  ie  OL 000001 


Target  Ground  Range  ( rve.  > 

Fl^jre  1-5.  emulative  Probdillity  of  Initial  Detection 
Detection  Criteria  le  Two  Coneeoutlve  BUpe 
Target*  -  Seerllng  Caee  I 

1  *q.e.  Radar  Croee  Section 
200  feat  Altitude 
480  knot  Closing  Rate 
Thrertnld  le  Fixed 

IVabdilllty  of  False  More  le  a  000001 
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Target  Ground  Range  (  a  m.  > 

figure  1-6.  Cumulative  Prcbcbility  of  Initial  Detection 
Detection  Criterion  is  Too  Coneecutive  Blips 
Targeti-  Generalized  Gaeeo  Distribution  <  K  «*  0.4  ) 
1  eg.  e.  Radar  Croee  Section 
200  feet  Altitude 
480  knots  Closing  Rate 
Threshold  is  Fitted 

Probcbiiity  of  False  Mare  is  0.000001 
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Target  Ground  Range  ( as.  ) 


Figure  1-5.  Cuevlative  Probability  of  Initial  Detection 
Detection  Criteria  ie  Two  Consecutive  Blips 
Tcrgeb-  Sterling  Case  I 

1  sq.s.  Radar  Cross  Section 
200  feet  Altitude 
480  knot  Closing  Rate 
Threshold  is  Fixed 

Probrfilllty  of  False  Alar*  is  0.000001 
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